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The bispectrum of scalar mode density perturbations is analysed for the strong regime of warm 
inflationary models. This analysis generalises previous results by allowing damping terms in the 
inflaton equation of motion that are dependent on temperature. A significant amount of non- 
gaussianity emerges with constant (or local) non-linearity parameter /jvl ~ 20, in addition to the 
terms with non-constant /jvl which are characteristic of warm inflation. 



_^ I. INTRODUCTION 

o 

Observations of the cosmic microwave background are consistent with the existence of gaussian, weakly scale 
I dependent, density perturbations as predicted by most inflationary models [IHl]. The amount of non-gaussianity 
, ^ produced by the simplest inflationary models is small and unlikely to be to be observable by the next generation of 
experiments, but this still leaves open the possibility that a slightly more exotic inflationary model could produce a 
measurable effect. 

One variation on inflation is the warm inflationary scenario [S] (see also [B] and the review [7|). The non-gaussianity 
,_i produced by a warm inflationary scenario with constant friction coefficient has been presented elsewhere [51. The 
bispectrum of the non-gaussianity is large, and it has a distinctive dependence on wave number. In this paper, we 
r \ shall examine the bispectrum of the non-gaussianity in the more likely situation where the friction coefficient is not 
constant. We shall see that the bispectrum has two terms, one which is like the previous warm inflation bispectrum 
and a new term which is typical of density perturbations with local non-linearities. 

Warm inflation is characterised by the rate of radiation production during the inflationary era. The radiation can 
Q affect both the homogeneous evolution of the inflaton field and the inhomogeneous fluctuations. If the radiation 
^ field has a strong damping effect on the inflaton dynamics, then we have what is known as the strong regime of 
c/3 warm inflation. If the damping effect is small, but the fluctuations are still influenced by radiation, then we have the 
I ^ I weak regime of warm inflation. In both strong and weak regimes of warm inflation, fluctuations in the radiation are 
transfered to the inflaton [5141 6| and become the primary source of density fluctuations. This is the most signiflcant 
CNJ difference between warm inflation and traditional cold inflation. 

^ The simplest warm inflationary scenarios assume that the radiation produced during the inflationary era thermalises 

at a rate faster than the expansion rate. This type of warm inflationary model is therefore rather restrictive. However, 
^2 the possibility of thermalisation occurring is enhanced by inflaton decay channels, which are naturally present in many 
supersymmetric theories, where the inflaton decays into light radiation flelds through heavy particle intermediaries 
[T7H2D] . In these models, the damping of the inflaton held is described by a friction coefficient F oc T'^, where T is 
the temperature. In particular, c ~ 3 at temperatures small compared to the heavy particle masses [211 HI]- The 
temperature dependence has been found to have a large effect on the size of the density fluctuations [33] , and we shall 
examine now how the temperature dependence affects the non-gaussianity. 

A feature of warm inflation is that a signiflcant amount of non-gaussianity is produced whilst the density fluctuations 
are still on sub-horizon scales. In the previous analysis, this was due to the non-linearity caused by the bulk velocity 
of the radiation. In the new analysis, there are terms in the non-gaussianity which are proportional to the parameter 
^ c. The only restriction on c for the consistency of warm inflation is that c < 4 [24] . This leads to effects which are 
large compared to slow-roll parameters. For comparison, the non-linearities produced by derivatives of the inflaton 
potential in warm inflationary models was looked at by Gupta et al. [251 126| . but their contribution to the non- 
gaussianity is governed by slow-roll parameters and it is tiny in comparison with the true non-gaussianity produced 
in warm inflation. 

Fluctuations in the cosmic microwave background provide an observational link to the density fluctuations at the 
surface of last scattering. We know, in principle, how to evolve these fluctuations from early times using, for example, 
the Bardeen variable C 127|. Observations can be compared to predictions for various moments of the probability 
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distribution of C- The most important of these is the primordial power spectrum of fluctuations (fc), defined by the 
stochastic average 

(C(ki)C(k2)) = (27r)3p^(fci)5'(ki + k2). (1) 

The bispectrum fc2, fcs), defined by 

(C(ki)C(k2)C(k3)) = {27rfB^{kuk2, h)S^k^ + k2 + ka), (2) 

can be used to examine the non-gaussianity in the density fluctuations. The normahsed amount of non-gaussianity 
in the bispectrum is described by a non-hnearity function fNL, deflned by 

where the 5/6 factor is convenient for cosmic microwave background comparisons [28j . Models with constant /^vl 
are often called local models because this type of non-gaussianity can arise from local non-linearities in the density 
perturbations. 

Non-linear evolution during the inflationary era can result in non-gaussianity appearing in the primordial density 
fluctuations. The amount of non-gaussianity produced by vacuum fluctuations in single-field inflationary models is 
typically around a few per cent PM5T] . and can be related to the standard set of inflationary slow-roll parameters 
[32] , This is small compared to non- inflationary effects. For example, the second order Sachs- Wolfe effect is expected 
to act as a source of non-gaussianity in the cosmic microwave background observations equivalent to /at l ~ 1 j33l I34j . 

Models of inflation with multiple scalar fields, acting as sources of density fiuctuations in the curvaton scenario 
[351 136) . or modifying the reheating phase of the universe '37], can produce a level of non-gaussianity above the 
foreground effects, even signiflcantly above the foreground for particular parameter choices. Non-gaussianity can 
also be produced by modifications to the kinetic part of the inflaton Lagrangian in D-brane models [38l |39] . This 
type of non-gaussianity is concentrated on equilateral wave- vector triangles fci, fc2 7 ^3, unlike the local form, which is 
concentrated on oblique triangles. 

The best observational limit on the non-gaussianity at present is from the WMAP seven-year data release [4^ , which 
gives —10 < /j^™' < 74 with 95% confldence for a constant (or local) component. The Planck satellite observations 
may have a sensitivity limit of around l/j^™'] ~ 5 — 10, depending on how well the signal can be separated from the 
galactic foreground 58] . Our prediction for /^£''' from the strong regime of warm inflation lies well above the Planck 
detection threshold in most models, but the most significant feature is the presence of a term in the bispectrum which 
could eventually provide a means to distinguish warm inflation from other sources of non-gaussianity. The prospects 
for observing this term have been discussed in Refs. f^T and f??. 

The paper is organised as follows. We begin in section II with a brief introduction to the notion of warm inflation. In 
section III, we introduce fluctuations of the inflaton fleld described by a Langevin equation and expand the fluctuations 
to second order. The limit of strong dissipation is introduced in section IV, and the bispectrum is calculated in section 
V. Some observational prospects are discussed in the conclusion. We use units with T% = c— 1. 



II. WARM INFLATION 



Warm inflation occurs when there is a significant amount of particle production during the infiationary era. We 
shall assume that the particle interactions are strong enough to produce a thermal gas of radiation with temperature 
T . In this case, warm inflation will occur when T is larger than the energy scale set by the expansion rate H . The 
production of radiation is associated with a damping effect on the inflaton, whose equation of motion becomes 

+ (3il + V)4> + = (4) 

where r(0, T) is a friction coefficient, H is the Hubble parameter and is the (j) derivative of the infiaton potential 

The effectiveness of warm inflation can be parameterised by a parameter r, defined by 



When r ^ 1 the warm inflation is described as being in the strong regime and when r ^ 1 the warm inflation is in 
the weak regime. 
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Consistent models of warm inflation [21] require a suppression of thermal corrections to the inflaton potential, so 
that the efi^ective potential separates into inflaton and radiation components 

F(0,T) = F(0) + p,(r), (6) 

where pr is the radiation density 

Pr = -^3*^^- C^) 

In this case, the time evolution is described by the equations 

+ (3i/ + r)0 + = 0, (8) 

Pr+^Bpr^V^^, (9) 

ZB^ = 4ttG (2V + 2pr + (10) 

During inflation we apply a slow-roll approximation and drop the highest derivative terms in the equations of motion, 

3Hil + r)<j) + V^ = 0, (11) 

4Hpr^r4,\ (12) 

37J2 = SnGV (13) 

The validity of the slow-roll approximation depends on the slow-roll parameters defined in |16j . 

1 ^vsV 1 ^vss\ „ 1 ^rsVA 



levrG \ V J ' 87rG[v )' ^ 8nG [ TV ) ^^^^ 

The slow-roll approximation holds when e<S^l + r,ri<^l + r and (5 <^ 1 + r. Any quantity of order e/(l + r) will be 
described as being first order in the slow-roll approximation. 

The temperature dependence of the friction coefficient T plays an important role in the present analysis. We 
parameterise this by a parameter c, 

TTt 

c=^, (15) 

where Tt denotes the T derivative of T. This parameter is not necessarily small, but a stability analysis of warm 
inflation shows that c < 4 for a consistent model [24l. 



III. FLUCTUATIONS 



Thermal fluctuations are the main source of density perturbations in warm inflation. Thermal noise is transfered to 
the inflaton field; mostly on small scales. As the comoving wavelength of a perturbation expands, the thermal effects 
decrease until the fluctuation amplitude freezes out jl3j . In the strong regime of warm inflation, this occurs when the 
wavelength of the fluctuation is still small in comparison with the size of the cosmological horizon. 



A. Inflaton fluctuations 



The behaviour of a scalar field interacting with radiation can be analysed using the Schwinger-Keldysh approach to 
non-equilibrium field theory [43l|44]. In flat spacetime, when the small-scale behaviour of the fields is averaged out, 
a simple picture emerges in which the field can be described by a stochastic system whose evolution is determined by 
a Langevin equation [IS]. This takes the form 

- V20(x, t) + T^{x, t) + v^ = (2rr)i/2^(x, t), (16) 

where is the flat spacetime Laplacian and ^ is a stochastic source. For a weakly interacting radiation gas the 
probability distribution of the source term can be approximated by a localised gaussian distribution with correlation 
function ^ ^ , 



{^{x, t)^{x', t')) = ,5(3) {x - x')5{t ~ t'). 



(17) 
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We shall restrict ourselves to this gaussian noise approximation. 

We can use the equivalence principle to adapt the flat spacetime Langevin equation to an expanding universe 
during a period of warm inflation by replacing ordinary derivatives with covariant derivatives in the cosmological 
metric with scale factor a and co- moving coordinates x". The Langevin equation will retain its local form as long as 
the microphysical and thermal scales in the problem are small compared to the cosmological scale [3117]. However, 
the rest frame of the fluid will have a non-zero 3— velocity with respect to the cosmological frame and we must include 
an advection term. Another alteration is that the coefficient of the noise term is changed slightly by the expansion 
(for details see ref ^J,to K = {2T^gTy^^, where Teff = T + H. 

The Langevin equation becomes [5] 

4){x, t) + m^ix, t) + TD(t) + V^- a-'^d^(l){x, t) = KC{x, t) (18) 

where is the Laplacian in the expanding frame and is the derivative along the radiation fluid. The correlation 
function for the noise, expressed in terms of the comoving cosmological coordinates, has the form 

{^{x,t)^{x\t')) ^a-^{2Tif5^'^\x~x')5{t~t'). (19) 

The fluid velocity components in a coordinate frame are u = (7,u"), where 7 is the Lorentz factor. With this choice 
of components, 

D(j) = -i^ + U"da(j3. (20) 

The inflaton will generate metric inhomogeneities, but with a suitable choice of gauge, these can be discarded on 
sub-horizon scales (see Sect. IIIC). We shall use a uniform expansion rate gauge. Eq. (18 1 applies on scales which 
are intermediate between the thermal averaging scale and the horizon scale. Later, we use a matching argument to 
extend the fluctuations to large scales. 

The analysis of the Langevin equation can be simplified by introducing a new time coordinate r = (aH)^^ and 
using the slow-roll approximation. We are led to the equation 

(j)"{x, t) - (37r + 2)r- V'(a;, r) - irr'^au'^dMx, r) - d'^(j){x, r) = Ki{x, r) (21) 

where a prime denotes a derivative with respect to r and we have kept only the leading terms in the slow-roll 
approximation. The noise term has been rescaled so that its correlation function is now 

{i{x, T)i(x', t')) = {2n)H^^\x - x')5{t - r'). (22) 

This equation is non-linear because F, it" and K depend on (j) and T . 

Now we treat the source term as a small perturbation and expand the inflaton field 

0(x, t) = 0(t) + (5i0(a;, t) -|- 52(j){x, t) + . . . (23) 

where 5i<f) is the linear response due to the source ^. Similarly, for the fluid velocity, 

u"(x,t) = <(x,T)-H<(a;,r) + ... (24) 

This expansion is substituted into the langevin equation. Only the zcroth order terms in the slow-roll approximation 
will be retained. 

The first two perturbation equations are 

5i(t)" - (3r + 2)T-^5i<t)' - d^5i(t> + 3i7- V^^^ d^r = Xi, (25) 
62(1)" - (3r + 2)t-i^20' - a2<520 + SH-\-^^S2r = SiK^ 

+ST-^6irSi(j)' - Sarr-^'^dJicI) - ir^H-^T-^621 (26) 

To leading order in the slow-roll approximation, the perturbations of r are determined entirely by the temperature 
dependence of the friction coefficient F. Since r cx and pr oc T'^, we obtain 

Oir = cr (27) 

^ Apr 2 \Apr ) ^ ' 
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Similarly, 

Before substituting these quantities back into the perturbation equations, it is convenient to replace perturbations by 
dimensionless parameters Cn a-nd e„, 

Cn = (30) 

On large scales, the parameters and ei become the Bardeen variables for an inflaton dominated and a radiation 
dominated universe respectively. Note that pr and </> are related by the slow-roll equation (|l2[). 



After substituting Eqs. (27) and (281 into the perturbation equations and converting to dimensionless variables we 
have 

C" - (3r + 2)t-^C - d^Ci + 3crT-^ei = (32) 
C^' - (3r + 2)t-\^ - d^C2 + 2,CTT-^e2 - SiKS, 

+ |c(4 - c)rT-'^ei + Scrr-^eiCi - 3orT"i<5aCi - irT-'^52l- (33) 

where K = 



B. Radiation fluctuations 



The dominant source of fluctuations in the radiation field is an inhomogeneous energy-momentum flux from the 
inflaton field. This transfer of momentum and energy into the radiation is described by an energy-momentum four- 
vector Qa, ElllH] 

Qa = -Tu''db(j)da(l). (34) 

where u is the 4— velocity of the radiation fluid, u= (7, u"). 

We shall model the radiation field by a perfect barotropic fluid with pressure p = wpr, and energy momentum 
tensor 

Tab = (1 + w)prUaUb + PQab, (35) 

and field equations 

Va'T"'' = Q". (36) 
The time and space components of the field equations are 

Dpr + {l+w)pr'^aU'' - Q, (37) 

wdipr + {l + w)prDu^ = Qi, (38) 
where _L denotes components perpendicular to u and the source terms are 

Q = -u'^Qa^TiDq^f, (39) 

Q„ = -r(z?(/))a„0. (40) 

As before, metric perturbations are small on sub-horizon scales and we can use the cosmological background metric 
with flat spacial sections. The fluid divergence is given by 

Wau" = + 3Hj + 7 (4f ) 

Indices are lowered with the background metric, so that Ua = a^u" and 7^ = 1 -|- UaU". 
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1, First order perturbations 



First order perturbations of Eqs. (37) and (38) give 



&lPr + 'iH{l + w)SiPr + [l + w)prdaU" = 5iQ, 
wdaSiPr + {I + Wjprilia + 3H{1 + w)prUia — SiQa- 



(42) 
(43) 



Perturbations in the energy and momentum fluxes ( 39 ) and ( 40 ) are caused by perturbations in the scalar field and 
perturbations in the friction coefficient, Eq. (27), 



SiQ = cHSiPr + 2T(f) Si(i 

SlQa = -T(j)da5l<l). 



(44) 
(45) 



From this point on we shall take a radiation fluid with w = 1/3. 

The dimensionless flux perturbation q„ — 5nQ/Q can be introduced, in addition to the perturbations Ci, and £i 
used previously. Using the time coordinate t = l/{aH), 



adau" = 3e'j — 12t + 3r ^gi, 
Ui — or Ui — ao Si — 3t . 

where 

91 = c£i - 2tC{ 
qia = -ardaCi- 

The velocity is given in terms of the inflaton and density fluctuations by 

ac)„< 3e'i - 6Ci - (12 - 3c)T-'^ei 

After eliminating the velocity, 

e" - (8 - c)T'^e[ + (20 - 5c)T-'^ei - ^d^si 
This equation was derived previously in Ref. |23l . 



(46) 
(47) 

(48) 
(49) 

(50) 



(51) 



2. Second order perturbations 



Second order perturbations of Eqs. ( |37[ ) and (38) give 

^2Pr + 3i/(l + w)62Pr + (1 + 'w)pr SqUj 

+DiSiPr + (1 + w)SiPr daU"^ + (1 + w)Hpr92 = S2Q (52) 
wdaS2Pr + (1 + w)prU2a + 3i7(l + w)prU2a 

+WUiaSiPr + (1 + w)SiPrUla + (1 + w)prDiUia — 52Qa " "la'^iQ (53) 

where Di = da and H92 — 3Hd2"/ + (527- In terms of dimensionless parameters and w — 1/3, 

adaU2 — 3s2 — 12x^-^62 + 3r^"^g2 

-4aei9„u? - SaDiSi - 3x^1527 + ^27' (54) 
uf - 5t~^u^ - ad°'e2 + 3x^1 = 

- 4£iuf + 8T-iei< + aDi< + St'^'^qi (55) 

The second order momentum flux is 

= -aTd"C2 - careid^Ci + ar^Cid^Ci, (56) 

and the second order energy flux is 

<72 = c€2 - 2rC2 - 5c(4 - c)e? - 2creiCi + r^Cf + 2arDiCi + 2(527. (57) 

An equation for the energy density fluctuations can be found by eliminating the velocity as before. We have done this 
using a computer algebra package, and the result is given in the summary section below. 
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C. Metric fluctuations 



We have been claiming up to now that the metric fluctuations play no role in the analysis. This can be checked 
using the gauge ready equations for metric perturbations from Ref. |48j . The new quantities we require at first order 
are the proper time perturbation ai, the scale factor perturbation the perturbed expansion of the hypersurface 
normal vectors ki and the perturbed shear of the hypersurface normal vectors xi- For convenience we pick a uniform 
expansion-rate gauge ki = 0, although the conclusions hold in for any reasonable gauge choice. 

When the metric perturbations are included and we take the leading order in the slow-roll approximation, the scalar 



field equation ( 32 ) and the fluid equation ( 47 ) become 

C" - (3r + 2)t-'^C' + fc^Ci + icTT-'^ei + r-^a'i - Br^^^i = i^l, (58) 
wf - 5t-i< = a9"(£i - ai) - ir'^ql (59) 

The Raychaudhuri equation gives the metric perturbation ai (see Eq. (14) in [JB] or Eq. (A22) in ^), 

£ 6r£i ^ AtC.[ + 6(1 + r)Ci 

(l + r)2fc2r2-3e/(l + r)+4e/(l + r)2 ^ 

Clearly ai « £i and ai << Ci, not only for sub-horizon scales fcr >> 1 but also for scales comparable to the horizon 
when /cr 1. The equations for the second order metric variation a2 follow a similar pattern, apart from the presence 
of a large number of extra quadratic terms in the first order perturbations. The effects of a2 are also suppressed by 
the slow-roll parameters. 



D. Summary 

We now have a full set of perturbation equations up to second order which can be used to determine the radiation and 
scalar field fluctuations £„ and Cn- The first order equations were solved numerically in Ref [23', and an encouraging 
feature of the numerical work was that the results agreed well with analytic approximations. We shall be making 
similar approximations to the second order equations in the next section. 

The first order equations are 

CI - (3r + 2)t-^C[ - d^Ci + icrr-^ei = (61) 
e'l - (8 - c)T-^e\ + (20 - bc)T-'^ei - \d'^£i 

-2C" + %T-\[ - a^Ci = (62) 

The fiuid velocity is irrotational at this order, and it can be determined by the equation, 

a(9„< = 3e'i - QCi - (12 - ?.c)T-^ei (63) 

The second order perturbation equations are 

C2 - (3r + 2)t-\'2 - 0^2 + 3crr-2£2 = S^ki 

+ |c(4 - c)rT-^el + Scrr-hiCi - SarT-^^daCi - ST-^rS2j, (64) 
£2 - (8 - c)T-h2 + (20 - 5c)T-^e2 - \d^£2 

-2C + 8r-iC^ - d\2 = j(^) +/^\ (65) 

where the source terms j^"-* are ordered so that increasing n implies decreasing size for large r, 

/I) = a2a"[(c£i-05aCi]-|fl'5"(3Ci+£i)9a(9Ci+2£i), (66) 
= (2c£i-2TC;)Cr + 27Cf -(24-2c)C;£'i + 6£f 
+32(3 - c)T-^Cex - (12 - c)(4 - c)r"^£'i£i 

+ i(17c-48)(c-4)T-2£2 ^aw^a^o; (67) 
= -11T-2J27, (68) 

where 

w = -7CJ + \e\ + |(6c- 35)r"i£i - llr^^Ci- (69) 
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IV. STRONG DISSIPATION AND SMALL SCALES 



The perturbation equations simplify considerably in the strong regime of warm inflation. In this regime, the thermal 
fluctuations are generated on small scales, compared to the cosmological horizon, and their amplitude approaches a 
simple power-law behaviour. This contrasts with vacuum fluctuations, which oscillate on sub-horizon scales and only 
freeze-out when they grow larger than the horizon. 

Fourrier transforms are defined with respect to the comoving coordinates x". 



C„(fc) = / (fxCnix)e 



:k-x 



(70) 



The fluid velocity decomposes into an irrotational scalar component u„ along k and a solenoidal component u^" 
perpendicular to k. Only the scalar component. 



(71) 



contributes to the scalar density perturbations at second order. In this section we shall we shall take r = r/3H >> I 
and concentrate on physical scales small compared to the horizon, where the wave number k >> aH. 



From Eq. ( 62 1 , the first-order density and scalar fluctuations are related by 

£i ~ -3Ci- 



This approximation can be substituted into Eq. (25 1. 



C" - (3r + 2)t^\[ + k\i - 9crT-\i = ki 



We solve this using a Green function. 



k 



G{kT, kT'){kT'y-^'' K^(t') dr', 



where the retarded Green function G(z, z') is given by 



G(z, z') = -z'^z'^ {J.+3ciz)Y,+3ciz') - J.+3c{z')Y,,+3,iz)) for z < z' . 



(72) 
(73) 

(74) 
(75) 



and v = 3r/2. A useful feature of the solution is that in the range 1 << kr << the fluctuation has an approximate 
power law behaviour, cx r~^^. This allows us to use 



We can use both the approximations ( 72 1 and ( 76 ) in the velocity equation ( 63 1 , to get 

kui « 36(I + c)r"iCi- 



At second order in the perturbation amplitude, the dominant terms for large kr in Eq. (65) are 

fc2£2 + 3fc2(2=j(i)(A:), 



where j^^-* is given by Eq. (66 1. The approximations (72) and (76) imply j*-^' w and 

£2 ~ -3C2- 



We use this together with the approximations ( 72 ) , ( 76 ) and ( 77 1 in Eq. ( 64 ) to get 

C^' - (3r + 2)r-\!, + k^^ - ^crr'^^ « 



where * denotes a convolution and 



As 
A, 



27 



c(4 -I- c)r. 



-108(1 + c)r, 
-?(l + c). 



(76) 
(77) 
(78) 
(79) 

(80) 

(81) 
(82) 
(83) 
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Note that, when the friction coefficient is temperature independent and c = 0, then only the and terms survive. 
These terms agree with the c = case investigated in Ref. [8], apart form a sign change in Ay caused by an error in 
the sign of the green function used in the earher work. 

The As term arises for the temperature dependence of the friction coefficient. This term is local in space, and we 
might therefore expect it to produce a local type of non-gaussianity. The Ay term arises from the fluid advection 
terms u^da in the perturbation equations, and this term is associated closely with the bulk behaviour of the radiation 
field. It distinguishes between vacuum and non-vacuum fluctuations. 



V. NON-GAUSSIANITY AND THE BISPECTRUM 



The small scale inflaton fluctuations develop a simple power-law growth well in advance of the time when they 
cross the horizon. During this stage the metric fluctuations are relatively small (in the uniform expansion-rate 
gauge). Eventually, the wavelength of the perturbations crosses the effective cosmological horizon and the wavelength 
grows to a regime where the metric perturbations become important. On large scales it becomes possible to use a 
small-spatial-gradient expansion (first formalised by Salopeck and Bond [JH] and later developed into the 'delta-iV' 
formalism [SUHS^] ). This approach allows us to define a perturbation ( which is is constant on large scales, even in 
the full non-linear theory. The bispectrum and the non-linearity of the density fluctuations can be approximated, to 
a reasonable accuracy, by matching the small and large scale approximations at horizon crossing. 

In a uniform spatial curvature gauge, the perturbation C is defined by 



c 



{x,t) 



(84) 



where H = H{(f>) and ip = (j>{(j)) are given by solving the slow-roll equations ( ll]|l3 ). When the inflaton perturbations 
are expanded as before in eq. ( 23 ) , we have 



(85) 



where (j) subscripts denote derivatives with respect to (j) and — H/4>- Note that, according to eqs. ( ll|l3 ), Cifxtt/C^, 
can be dropped because it is first order in the slow-roll expansion. We therefore have 



C = Ci + C2 + 



(86) 



where the Cn are identical to the parameters defined in Eq. ( 30 ) of the previous section. 

The power spectrum and bispectrum of the density perturbations were defined in the introduction. The first order 
perturbations Ci are gaussian fields and their bispectrum vanishes. The leading order contribution to the bispectrum 
must therefore include a contribution to C from the second order perturbation, 



^ (Ci(ki,r)Ci(k2,T)C2(k3,T)) « (27r)3i?^(A;i,fc2,fc3)'53(ki +k2-Hk3), 

cyclic 



(87) 



where 'cyclic' denotes cyclic permutations of {ki,k2,k3}. The second order perturbation can be obtained by solving 
Eq. using Eq. ([74|. 

We shall split the second order perturbation C2, into three parts C2 = Cs + Cv + Cr, where each part is sourced by 
the corresponding terms with coefficients Ag , Ay or A^ in Eq. ( 80 ) . Beginning with the Cs term, we have 

C - (3r + 2)t-iC + k\s - 9crT-\s = AsT'^Ci * Ci- 



The solution can be obtained using the green function (75), 



/oo 
kdr'GikT, fcr')(fcr')-i-''^Ci*Ci 



(89) 



The contribution to the density fluctuation bispectrum from is 
{2TrfBg{ki,k2, fc3)'^'(ki + k2 + k3) = 
As /"MT'G(fc3T,fc3T')(fc3T')-i-'''(Ci(ki,r)Ci(k2,T)Ci*Ci(k3,T')) 



(90) 



cyclic 
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Now we can use Eq. ( 74 1 for Ci and decompose the stochastic average of the noise terms into products of the correlation 
resuh is that 

T>4: poo 

B,^2As -3-3 / fc3dr'G(fc3T,A:3r')(fc3T')-^-'''F(A:ir,fcir')i^(fc2T,fer'), (91) 



function ( 22 ) . The resuh is that 

-4 



cyclic 

where 



/"OO 

F{kTl,kT2) = fc / dT'G{kTl,kT')G{kT2,kT'){kT') 

J T-2 



(92) 



This integral can be evaluated analytically when v is large. The integral is expressed in terms the power spectrum 
P(^{k,T) in appendix [A} and gives 

fc3dr'G(fc3T, ksT'jik^TT^-^"'^-. (93) 

cyclic 

This integral can also be found in appendix |Xj leaving 

2A 

Bs = ^/(C) ^ P^{k^,T)P^(k2,T), (94) 
cyclic 

where we have used u = T/2H. The factor /(c) is, 

/(c) = 1 - (2r/i/)-3-/2. (95) 



We can usually take /(c) « 1 unless c is very small. The value of As was found in the previous section (81 ), so that 

= 3(4 + c)/(c) ^ Pc{k,,T)P<:ik2,T). (96) 

cyclic 

This part of the bispectrum has a local form, as expected, and the value of the non-linearity parameter for Bs alone 
is, 

/^T' = ^(4 + c)/(c). (97) 

Remarkably, this is only weakly dependent on F. However, some caution needs to be exercised due to the rather crude 
nature of the matching procedure applied at the horizon scale, where the approximations used in Sect. IV break 
down. 

In Ref. 1231 , an alternative approximation scheme was introduced to cover fluctuations at the horizon scale. This 
approximation implied that the density fluctuation power spectra approach constant values at the horizon scale, with 
P^(fci,r) = Pc(fci) oc k~^ (or at first order in the slow-roll parameters). In the present context, this suggests 

that the small scale approximation will extend to the horizon scale when the bispectrum is expressed in terms of the 
power spectrum, as we have done above. 



The bispectrum reaches its largest values for oblique wave- vector triangles, with ki << fc2 ~ fcs (up to permutations 
of the sides). However, if the triangle becomes too oblique, then the result no longer applies because it is not possible to 
have fciT and k2T in the same range where the approximations are valid. The solution [5] is to cut off the bispectrum 
when k2/ki > r = T/3H (or similarly for any other pairs of sides). There is some residual dependence of the 
bispectrum on F due to this effect. The bispectrum Bs with a cut-off for ^2/^1 > r is plotted in figure [l] 

The other contributions to the bispectrum are obtained by following the same steps. The result for By is 

By = ~{l + c)/(c) ^ (fcr' + kf)ki ■ k2 Pc(ki)P^{k2), (98) 

cyclic 

Note that, in the limit c — > 0, the function /(c) ^ (3c/2) ln(F/2i/) and the result agrees with Ref. [5], apart from 
the overall sign mentioned earlier. The bispectrum is concentrated on oblique triangle shapes as before, but this 
time the contribution to /nl is strongly dependent of the shape of the wave vector triangle and even changes sign 
as the triangle becomes more oblique. The dependence of the bispectrum By on ki and ^2 for r = 20 is plotted in 
Fig. [2j The contribution from the squeezed triangles is reduced relative to the local case, but the squeezed triangle 
contribution rises for larger values of r 
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FIG. 1: A contour plot of the truncated local bispectrum Bs as a function of ki and ^2. The plot shows 
Bs{xi, X2, l)x\x\ tanh^(ra;i) tanh^(ra;2)/-Bs(l, 1, 1), where xi = ki/ki, X2 = fc2/fc3 and r — 20. 




FIG. 2: A contour plot of the truncated local bispectrum Bs as a function of fci and ^2. The plot shows 
Bv(xi, X2, l)x\x2 tanh^(ra;i) tanh^(ra::2)/i3i,(l, 1, 1), where xi = ki/k^, X2 — k2/k-i and r — 20. 



The final contribution to the bispectrum depends on the integral (A18) 
factors of r/3-ff. It is only relevant when c ' 
bispectrum is effectively the sum of two terms 



This contribution is suppressed by 
0, and even then this contribution is dominated by B^,. The total 



B 



6 



flocal 

Jnl 



cyclic 



6 



'adv 
NL 



cyclic 



fc2"')kl •k2Pc(fci)Pc(fc2), 



(99) 
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where f'iif!' is named after fluid advection terms 



fadv 

Jnl ■ 



The negative sign in front of /^'^ is chosen so that, for equilateral 



NL 

triangles, /atl 

The values of /j^£"' and are plotted in Fig. Sl The contribution to the bispectrum from By dominates at 
small values of c, but then decreases, whilst the localcomponent Bs becomes larger as c increases. The momentum 
dependence makes the contribution to B from By harder to detect, but since this component is characteristic of warm 
inflation it would be desirable to try to isolate this component in future CMB experiments. 



60 p— 
50 - \ ■ 
40 - 
30 - 



jlocil 
^adv 




FIG. 3: A plot of fNL for equilateral triangles against c. The curves represent two different contributions with different 
dependence on wave numbers. In this plot, T/SH — 10 and F oc T*^. The contribution from /at'l' should be multiplied by the 
values in Fig. 1 for other triangle shapes, and experiments which search for constant (or local) fj^i^' will see predominantly 

/Local 
NL ■ 



VI. CONCLUSIONS AND OUTLOOK 



The main results of this paper are given in Fig. [3] which shows the amplitudes of two contributions to the bispectrum 
for primordial fluctuations which originate in the strong regime of warm inflation, with friction coefRcient T (x T'^. 
Contemporary models of warm inflation typically have c = 3 [20] . The local component with amplitude /^£°' is easier 
to detect, and would most likely be seen flrst. 

It would be very desirable to search for the component of the bispectrum with amplitude /^"^ , since this component 
is a characteristic feature of fluctuations which originate from a non-vacuum source. Some idea of the difficulty of 
measuring both components of the bispectrum can be gauged by considering an ideal experiment where the only 
source of noise is the cosmic variance. We can parameterise the bispectrum as we did in Eq. (99 1, 



b = Y,Ub\ 



(100) 



where i — 1,2 for /)^£° of fff-^ respectively. The observations give a set of spherical harmonic components Bf^'f^i^, 
with an upper limit / < Zmax- The estimator for the parameter /j is 



dCi^Ci^Ci^ 



(101) 



j hhh 

where Ci are obtained from the power spectrum, and the Fisher matrix 



^hhh^hhh 

QCi^Ci^Ci^ 



(102) 
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The Fisher matrix for the relevant models was evaluated in Ref. [H] , where it was found that the standard deviation 
of the estimator is around 5 times larger than the standard deviation in the estimator /^£"'. For Planck, the 
detection limit for /j^™' is expected to be around 5 — 10, depending on how successfully the backgrounds can be 
removed. This would imply that Planck would only be able to detect the presence of fffl^ if the value was at least 
25. The problems arising from cosmic variance could be overcome by taking higher resolution surveys, but then the 
issue of backgrounds contributing to /^"^ has to be addressed. 

There is another regime of warm inflation, the weak regime, where the approximations we have used here all break 
down. Nevertheless, the dependence of the bispectrum on the wave-vector triangle will contain the same components 
that we found here, with the possibility of two extra components appearing from terms which we have been able to 
discard in the strong regime of warm inflation. In the weak regime we have to resort to numerical investigation of the 
second order perturbation equations to determine the bispectrum fully, and this is something we hope to report on 
at a later date. 



Appendix A: Integrals 

We begin with an approximation to the integral 

/"OO 

I{kT)^k dr' G{kT,kT'){kT')-^-^''-^'' (Al) 



where v ~ T/2H and the retarded green function G is given in Eq. (75 1. The leading contributions to the integral 
for large and fixed t come from 

/-oo 

/(fcr) « --(fcT)-3<=(fcr)'^V^,(fcr) / J^,{z) z'^-'^'-^'^dz, (A2) 
where i/' = ly + 3c. We also have 



kr 



which give 



YAz)--z-'-TnM, Mz)^z''-^-—. (A3) 
TT TR[iy + 1) 



bcv 



The approximations used above fail when c is small. In that case, there is a saddle point in the integrand which comes 
to dominate the result. The easiest way to obtain the contribution from this saddle point is to obtain the result for 
c < and extend to c > 0. When c < 0, 

poo 

I{kT)K--{kT)-^''{kTY'Y^,{kT) J^,{z)z-^-'''-^''dz. (A5) 
2 Jo 

The integral, 

J,,(z)z dz^I rfl(^' + l + 3c/2)- ^^^> 



For large v and small c. 



I{kT)^-^{kT)-^^{Av)-^^l\ (A7) 
bcv 



Combining the two asymptotic results together gives 

I{kr) ^ ^{kry (l - {4iy)-'''/') • (A8) 

The next integral is 



/oo 
dr'GikT, kT')G{kT, kT').{kT'f~'^'' 



(A9) 



This integral arises when we calculate the power spectrum using Eq. ( 74 1 , 
The leading contributions to the integral for large v and fixed r come from 



4 Jo 
where v' = v -\- 3c. This gives a standard Schafheitlin integral, 

1 r;^(i.~i)rji(3c + 3/2) 



J^,{zYz^~"'dz 



Tr{v - 1/2)Tr{2v + 3c - 1/2) 



where T n is the gamma function. Hence, 

F{kT) r- 



TT r/j(3c + 3/2) / 2v 



3c 



2^2 



32j/ r7j(3/2) V^^^^ 

for fcr << i^^/^. The 0(i/^^) terms come from the expansion of Yi,{kT). 
The next integral is 



F{kTi, kT2) = k / dT'G{kTi,kT')G{kT2, kT'){kT') 



This is evaluated in exactly the same way as F{kT), 

F{kTl,kT2) 



3c 



F{kn). 



We can also use Eq. (AlOl to express this in terms of the power spectrum, 

k^FikTukT2) - 



3c 



Another integral of this type is 

F(fciT, k2T) - (fcifca)^/' / dT'G{k,T, k^T')G{k2T, A^st') (^l t') l"'" (^2 



This time, for large v, 



F(klT, k2T) 



2kiki 



(3c+3)/2 



{2v)-^^/^F{k,TMT). 



This integral is negligible compared to F{kiT, k2T) 
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